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sponding redshift of light received from them can only be explained by an expansion 
of the space between them and us. In this paper, for the presently favored case of a 
universe without spatial curvature, it is shown that this interpretation is restricted 
to comoving coordinates. It is proven by construction that within the framework 
of general relativity other coordinates exist in relation to which these phenomena 
can be explained by a motion of the cosmic substrate across space, caused by an 
explosion like big bang or by inflation preceding an almost big bang. At the place 
of an observer, this motion occurs without any spatial expansion. It is shown that 
in these "explosion coordinates" the usual redshift comes about by a Doppler shift 
and a subsequent gravitational shift. Making use of this interpretation, it can easily 
be understood why in comoving coordinates light rays of short spatial extensions 
expand and thus constitute an exemption from the rule that small objects up to the 
size of the solar system or even galaxies do not participate in the expansion of the 
universe. It is also discussed how the two interpretations can be reconciled with each 
other. 
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INTRODUCTION 



After its discovery, the redshift of light from far away galaxies was first explained in the 
context of special relativity and attributed to a Doppler shift caused by an outward flight 
of the galaxies in a preexisting invariable space. It was first observed by V. M. Slipher in 
the years since 1912 (see e.g. and in 1918 C. Wirtz interpreted it as being due to a 
general recessive motion of galaxies {2]. Use of the Doppler formula led to good agreement 
with observational results and revealed the validity of the Hubble law that, in fact, was first 
derived in 1927 from general relativity by G. Lemaitre Only in 1929 it was formulated 
by E. Hubble in the context of astronomical observations [4]. Hubble's interpretation of 
the cosmological redshift as a Doppler shift became the generally accepted view, and in 
accordance with it the big bang was considered as a giant explosion. 

The successful general relativistic formulation of the basic cosmological equations in co- 
moving coordinates, first by A. Friedman [fj] in 1922, later independently by Lemaitre jsj 
in 1927, and 1936 supplemented by a rigorous derivation of the corresponding metric by 
H. P. Robertson [6] and A. G. Walker [7], brought a second interpretation into play. Ac- 
cording to this the recession of galaxies from a distant observer is not caused by a motion 
relative to their spatial environment but rather by an expansion of the space between them 
and the observer. Correspondingly inflation and the big bang are sometimes denoted as 
space explosions. 

For a long time the two interpretations coexisted because in those days the most distant 

observable galaxies gave rise to very small redshifts z only, and at small z the velocity - 

redshift relation is the same for both interpretations. However, at least since 1998 in the 

f] 

context of the Supernova Cosmology Project [8j more distant cosmic objects (supernovae) 
were observed and discrepancies between the two interpretations became evident, the em- 
phasis has completely shifted to spatial expansion. By now even for small z the recession 
of galaxies is exclusively attributed to spatial expansion whereas the explosion perspective 
is completely ruled out. This has become the standard doctrine in research articles (see 
e.g. {9] and jlO|). in textbooks (see e.g. and 16]), and in Wikipedia 24], the inter- 

pretation as an explosion sometimes being denoted as a popular or common misconception 
(see e.g. jsj], []i|-[l5] or 16], page 28). In the world wide web a multitude of contributions 
arguing against the interpretation as an explosion can be found (see e.g. |17l]-j23|). 
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In 1934 W. H. McCrea and E. A. Milne 



25] showed that the Friedmann-Lemaitre equa- 



tions describing the dynamics of the cosmic scale factor a(t) in general relativity follow in 
exactly the same form from Newton's laws of motion and of gravity (more precisely: from 
the corresponding fluid equations by Euler). In this so-called Newtonian cosmology the 
galaxies are located at the positions r(t) = [a(t)/a(to)] r(t ) where r is the radius vector in a 
Euclidean space and r = is our position in the universe, and they move with the velocities 
v = H(t) f (where H (t) = a(t) / a(t)) relative to a preexisting and invariable space. Near the 
origin the velocities of galaxies and the gravitational field produced by them are so small 
that the Newtonian equations of motion should asymptotically coincide with equations ob- 
tainable from general relativity. This suggests that also in general relativity there should 
exist coordinates in which the galaxies are moving across space. 

In this paper, for the presently favored case of a universe without spatial curvature, it will 
be shown that the seemingly contradictory interpretations, motion relative to an invariable 
space and recession due to expansion of space, are both possible and not in contradiction 
with each other. It will become apparent that they are relative views restricted to special 
sets of coordinate systems which are related to each other by simple transformations. 

It is clear that in the transition from one coordinate system to another also the inter- 
pretation of physical phenomena can change. Since in general relativity a huge variety of 
different coordinate systems is available, it may at first glance appear almost trivial to find a 
coordinate system which allows for a prescribed interpretation - motion across space instead 
of expansion of space in our case. There are, however, several restrictions making things less 
trivial. In order to remain within the framework of general relativity, the signature of the 
metric of all admissible coordinate systems must be the sameQ In addition, conditions for 
excluding radial expansion must be imposed. This will lead us to a set of differential equa- 
tions to be solved under the observance of specific boundary conditions. It is not obvious 
that solutions of this problems exist, and even when they exist, it is still not granted that 
the corresponding metric has the appropriate signature. A simple example may demonstrate 
the problem. 

In Newtonian mechanics one can go from the Cartesian coordinates of an inertial system 
to any other coordinates and still remain within the framework of the theory, provided 
1 In this paper SI units are used whence the speed of light is c. Furthermore the signature of the metric is 

(+, " " ")■ 
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that pseudoforces are admitted. Specifically one can go to a rigidly rotating coordinate 
system. The fact that beyond a certain radius this system rotates with superluminal velocity 
poses no problem in Newtonian mechanics. In general relativity, however, rigidly rotating 
coordinates are not admissible because in the superluminal region the metric assumes the 
wrong signature. A second example is presented in footnote [2] 

The main body of this paper consists of proving the existence of explosion coordinates. 
This is achieved in a constructive way deriving them by means of a transformation from 
Robert son- Walker coordinates. In section |T] the determining differential equations and 
boundary conditions are specified, formal solutions are derived and properties of them as 
well as conclusions following from them are discussed. In Section HT1 explicit solutions for the 
most important special cases of cosmic evolution are deduced, thus bringing to an end the 
envisaged proof of existence. As already mentioned the interpretation of physical phenom- 
ena like the redshift of light from far away galaxies must be adjusted. It will be shown that 
the latter can exactly be described by a Doppler shift and a subsequent gravitational shift. 
In this context it will be discussed why and how a spontaneous Doppler shift at the place 
of emission can be reconciled with the fact that in co moving coordinates "... the increase of 
wavelength from emission to absorption of light does not depend on the rate of change of 
[the cosmic scale factor] a(t) at the times of emission or absorption, but on the increase of 
a(t) in the whole period from emission to absorption" 12]. In section [Til it is also discussed 
to what extent the cosmic explosion provided by inflation or a big bang can be compared 
with usual explosions. Section II III deals with some applications and extensions. 

It is clear that the purpose of this paper cannot be the replacement of the usual ap- 
proach with Robertson Walker coordinates and the corresponding interpretation (also, see 
Section HV|) . especially since it starts off with solutions obtained in them. Rather the paper 
is meant to provide a supplementation which might even turn out to be useful in specific 
cases. 



I. SYSTEMS WITHOUT RADIAL EXPANSION GENERAL THEORY 

In the case of a spatially uncurved universe the square of the line element in (comoving) 
Robert son- Walker coordinates is 

ds 2 = c 2 dt 2 - a 2 (t) (dr 2 + r 2 dtt) with dtt = d$ 2 + sin 2 tf dtp 2 . (1) 
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(r is dimensionless and a(t) has the dimension of a length.) The radial expansion of the 
universe is expressed by the time dependence of length elements, e.g. dl r = a(t) dr in radial 
direction. The underlying coordinates t, r, •& and <p are called expansion coordinates in this 
paper. We are looking for a transformation t, r, tp — > r, p, ip to new coordinates r, p, d 
and (p, the explosion coordinates, for which the square of the line element is given by 

ds 2 = c 2 g 00 (p, t) dr 2 - dp 2 + g n (p, r) dtt . (2) 

(In contrast to r, p has the dimension of length.) The radial length element dl p = dp is time 
independent whence in explosion coordinates there is no radial expansion. As we shall see 
instead of this there is a radial motion. Since $ and ip remain unchanged, we have 

t = t(p,r), r = r(p,r). (3) 

Other than in Newtonian cosmology the space is not completely invariable but may involve 
angular expansion because go, can be time dependent. A corresponding expansion of volumes 
will be discussed at the end of this section. 

1. Derivation of explosion coordinates and corresponding metric 
From eqs. we get 

dt = t p dp + t T dr , dr = r p dp + r T dr 

where t p denotes the partial derivative of the function t(p,r) with respect to p etc.. With 
this, from eq. ([I]) we obtain 

ds 2 = (c 2 t 2 p - a 2 r 2 p ) dp 2 + (c 2 t 2 T - a 2 r 2 T ) dr 2 + 2(c 2 t p t T - a 2 r p r T ) dr 2 + . . . . 

In order for this to assume the form of eq. (j2J), the equations 

2/ / 2 2j.2 22 i j.2 a T t 2 2 ( a\ 

c t p t T = a r p r T , ct p -a r p = -1, g 00 = t T — , g n = -a r (4) 

must be fulfilled, or, equivalently 

Jl + cHl c 2 t t 

P <t) ' Tt ± a{t )JlT^ p - [) 
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In order to make the results in expansion and explosion coordinates comparable we impose 
the following conditions: 1. the origins of the two coordinate systems permanently coincide, 
and 2. the coordinate times at the origins are the same, 

1. p = for r = , 2. i = r at p = . (6) 

Eqs. (jHJ) are two equations only for the determination of the four derivatives r p , r T , t p and 
t T , since as soon as the transformation functions ([3]) have been determined from eqs. (j3J), 
the last two of the eqs. P| only serve for the evaluation of g o an d gn- For eqs. ((Sj) to have 
solutions, the integrability condition r pr = r rp must be satisfied. The evaluation of it results 
in a nonlinear second order differential equation for the function t(p, r), 

a(t) t pp + (c 2 t 2 p + 1) = . (7) 

Once a solution t(p, r) is found, r(p, r) can be determined from eqs. (jSJ). After multiplication 
with 2a t p eq. d7j) becomes 



2a(t) d(t) d 



2a 2 (t) t p t pp + 2a(t) d(t) t 3 p + tp = __ a 2 (t)(t 2 p + 1/c 2 ) 

This equation can be integrated once to yield 



0. 



t p = ±-jF{r)/a\t)-l (8) 



(9) 

r=const 



and once more to yield 

P = G(T)±c[ t dt ' 

J o ^F{T)/a 2 {t') - 1 

where F(t) and G(r) are integration functions. (Note that the algebraic sign in eqs. (jHJ) 
and (jSJ) can be chosen independently from that in the eqs. (jSJ).) The latter can be chosen in 
such a way that additional conditions are satisfied. Using the second of the conditions (jSJ), 
from eq. (jUJ) we get 



G(r) = T c T 
Jo 



(10) 

T=const 



'o ^F{t)/o?{V) - 

Inserting eq. (|H]) and the second of the eqs. (jSJ) in the third of the eqs. (@J yields 



Since under observance of both conditions © the origins of the two coordinate systems 
permanently coincide and since t — r there, also the metric times must be the same there. 
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(Although this statement needs no separate proof, it will later be validated as a test in 
special cases.) As a result of this from eqs. (JT]) and (j2J) we get 



9ao(p = 0,t) = 1 



(12) 



With this and t T = l for p = since t = r for p = 0, from eq. ( JIT]) we finally obtain 



F(t) = a 2 (r) . 



(13) 



Here, a(r) stands for the function a(t) with t replaced by r. With this and eq. (fTUj) eq. 
becomes 



p = ± c 



a 2 (r)/a 2 (^ 



(14) 



r=const 



This equation implicitly determines the function t = t(p,r). Inserting eqs. (fT3|) and (jSj) in 
the first of the eqs. (jSJ) yields 



"a 2 (t) 



and r = i?(r) ± a(r) / 

J o 



dp' 



o a 2 (t(p',r)) 



(15) 



T=const 



where -R(t) is an arbitrary integration function. From the first of the conditions ([6]) it follows 
that R(t) = 0, so finally we get 



r(p,r) 



air 



dp' 



o a 2 (t(p',r)) 



(16) 



T=const 



(The minus sign could be excluded since we must have r > 0.) 

The eqs. (Tlli) and ( fl6i) determine the solutions of the eqs. (00) for the transformation 
functions r(p, r) and £(p, r) only implicitly and require the knowledge of a solution a(t) 
of the Friedman-Lemaitre equations. Since according to eq. fHM their range of validity is 
restricted to a(r) > a(t), it is at this point not yet clear whether for each a(t) also explicit 
solutions exist and in addition lead to the right signature of the metric. Therefore, the proof 
of existence of solutions is not yet complete. Their determination will in general require 
numerical methods. However, in special cases also analytical solutions can be obtained. In 
section [III the full solution for two relatively simple but nevertheless representative cases will 
be derived. 
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2. Flight velocity of galaxies 

In explosion coordinates the radial position of a galaxy (or an element of the cosmic 
substrate) is given by r(p, r) = const, whence p = p(r) and 



r T ca 2 (t)t T y/a*{T)/a*{t)-l 



(For the last step eqs. (J5]), ([8]) and ( TT3|) were used.) This means that the galaxies are moving 
across space. p(r) is the coordinate velocity, and from it the physical velocity v is obtained 
according to ^ dp dp p(r) o(r) 



y/goodr y/g^dr a(t)t T 
where at last eqs. (TTTT) and (|T3l) were used. Insertion of eq. (|T7|) finally yields 

v = c^l-a 2 (t)/a 2 (T). (19) 

Thereby a freely disposable sign was chosen such that v is positive. 



3. Redshift of light from far away galaxies 

1. In expansion coordinates, the redshift z of light emitted from a galaxy at the radius 
r em at the time t em and received at the origin r = at the time t is given by 

^ + 1 = 3r = sr\ < 20) 

where X em or Ao is the wavelength of the light at emission or reception respectively. The 
radial propagation of light is described by cdt = a(t) dr whence 

ft° dt' . . 

r em = cl W) . (21) 



With this equation the time t of reception can be expressed in terms of the time t em and 
the location r em of emission. 

2. In explosion coordinates for an observer fixed at the distance p = p e m from the origin 
the light emitted from a galaxy flying past him at the radial velocity v undergoes the 
longitudinal Doppler shift 

\ p fl+v/c\ 1/2 l + v/c 



A, 



(22) 

1-v/cJ ^Jl-v 2 /c 2 
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^ (l + y/l-a?(t)/a?(T)") , (23) 



The physical velocity of eq. (fl9|) is the one that must be inserted in the Doppler formula, 
and because unlike the recessional velocity a(t) r in expansion coordinates it cannot exceed 
the speed of light, it is always in the range of validity of this formula. With it eq. (1221) 
becomes \ 

A e m Ct[t) 

where r = r em or t = t(p, r) = t em are the times of emission in explosion or expansion coordi- 
nates respectively. On its way from p = p em to the observer at p = the Doppler-shifted light 
experiences an additional shift in the gravitational field acting in the system of explosion 
coordinates. For simplicity we assume that this field is time independent as is true for the 
special cases to be studied below. Then, according to general relativity, the ratio of the 
wavelength Ao observed at p = and r = t to the wavelength X p of the Doppler-shifted 
light at p em — > p and r em — > r is 

Ao = 9oo{p = 0,r = r ) 
A P \ 9oo(p,r) 

Making use of the eqs. (fTTj) -( ri"3"l) . from this we get 

Aq 1 air 



(25) 



X p ^oo(p,r) <t)tr' 
Combining eqs. (|23|) and (|25|) we finally obtain 

Since observers at the origin use the same metric time in both coordinate systems, they 
also measure the same frequency or redshift of light that is emitted by far away galaxies. 
Therefore the result ( 1261) must be the same as the result provided by eqs. ()20l)-(l2TT). and 
although not evident this coincidence does not need proof. Nevertheless it is illuminating 
to see it verified in the special cases considered in section [Til 



4- Volume expansion rate 

One could be tempted to assume that the elimination of radial expansion in explosion 
coordinates leads to increased angular expansion. Indeed the structure of gn(p,r) in ex- 
plosion coordinates (see eq. (1321) for example) suggests that there could still be a volume 
expansion due to a time dependence of azimuthal distances. In order to clarify this issue 
we calculate an expansion rate that involves both radial and angular properties, the volume 
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expansion rate E = (dAV/dT)/AV. In this AV^ is the volume of a spherical shell with 
infinitesimally small thickness Ap, and T = ^Jgm r is the time measured on clocks in the 
system of explosion coordinates whence 

A t r rr. — j_ \ / 



AV y/goodr 
According to eq. (j2J) and the last of the eqs. (jlj) 

AV = An g a Ap = An a 2 (t) r 2 Ap 

with t = t(p, t) and r = r(p,r). From this, for fixed p we obtain 

1 dAV _ 2[a(t)rt T + a(t)r T ] 
AV dr a(t)r ' 

and using eqs. (TTTT) and (|T3l) we finally get 



E= w(m + IL y (28) 

a(t) \a(t) rt T j y ' 

For obtaining the expansion rate E at p = we consider the volume 

V = An f a 2 {i)r 2 dp' 
Jo 

of a sphere of radius p. Expanding around p = we have 

r(p, r) = r(0, r) + r p (0, r) p + . . . = r p (0, r) p + . . . 

because r(0,r) = according to the first of the conditions (jSJ). Furthermore, according to 
the first of the eqs. (|T5l) we have 

1 



r P (0,T)- 



p=o a l r J 



a 2 (t) 

since £ = r for p = according to the second of the conditions ([6]). In consequence 



r(p,r) = -^- + 0(p 2 ). 
a{T) 

Finally, we get 

a(t(p, r)) = a(t(0, r)) + O(p) = a(r) + 0(p) 

Altogether we have 

, An p 3 



V = An [ P [p> 2 + 0(p' 3 )] dp' = 4^ + 0(P 4 ) • 
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From this we get dV/dr\ p= o = and E\ p=0 = 0, in explosion coordinates at the origin there 
is no volume expansion. Since in an expanding universe each point can be chosen as the 
origin, in explosion coordinates at each point for a local observer there is no expansion. The 
expansion around points at some distance from the origin, described by equation ( )28|) . can 
therefore be regarded as virtual. 



II. SYSTEMS WITHOUT RADIAL EXPANSION SPECIAL CASES 

For simplicity our analysis of special cases will be restricted to situations in which the 
gravitational field is static in explosion coordinates. Since Newtonian cosmology is a very 
good approximation in the neighborhood of the origin r = 0, the corresponding conditions 
can be determined with it. According to the Hubble law v = H(t)r the accelerating (or 
decelerating) field, equaling the gravitational field (an inflation field being included), is 

v(t) = H{t) r + Hr[t) = (H(t) + H 2 )r . 

It becomes time independent for 

Hit) +H 2 = 

In the first case, we obtain 





A 2 = const 
-A 2 = const . 



„, s ait) 1 . s 

H(t) = -YT = - aU) = at 

where a = const, i.e., the universe is expanding with constant velocity. In the second case 
two solutions Hit) exist, firstly H = A, i.e. a universe with constant expansion rate, and 
secondly 

H(t) = A coth(Ai) ->■ a(t) = a sinh(At) . 
In the third case we obtain 

H(t) — A tanh(At) ->■ a(t) = a cosh(At) . 



In the following we shall only consider the cases of constant expansion velocity and constant 
expansion rate. 
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A. Universe with constant expansion velocity 

For a long time the recent evolution of the universe up to the present state was considered 
to be best represented by a model that assumes the domination of matter and exhibits a 
decelerated expansion with a(t) ~ t 3 ^ 2 . Only shortly before the end of the last century it 
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. The case of a universe 



was detected that the expansion is rather slightly accelerated 
without acceleration or deceleration, a(t) = at with constant a, lies in between, represents 
a fairly good approximation - which is best at the time of the transition from decelerated 
to accelerated expansion - and enables a full analytic solution of our problem. 



1. Derivation of explosion coordinates 
With a(t) =at the integral in eq. ( fl~4l) can readily be evaluated and is given by 

Jr y/a?(T)/a?(t') - 1 2 J - Vr 2 - t' 2 
With this, resolving eq. ( |T4l) with respect to t and choosing a sign such that t > 0, we obtain 

1 



t = - ^c 2 r 2 - p 2 . (29) 
With a(t) = at and this result eq. ( 1T6I) becomes 

r = ar (" dp ' = — [" dp ' = — In ° T + 9 (30) 
Jo a 2 t' 2 (p',T) a Jo c 2 t 2 — p' 2 2a ct — p 

The transformation t, r, if — > T,p,"d,ip, providing a line element of the form (T2]), is given 
by eqs. ( 129]) and ( 130]) . and the inverse transformation is 

. ar . ar 

r = t cosh — , p = ct smh — . (31) 

c c 

Inserting eqs. (1291 and (1301 in the third and fourth of the eqs. (jl]) yields 

*>=!, 9 n = -^^(ln^y. (32) 

4 y ct — p J 

(In the derivation of the first equation, for the sake of brevity eq. ( fl2i) was used. One can 
easily make sure that the same result is obtained without it.) Inserting the results (132]) in 
eq. ([2]) yields 



,2_ 2,2 ^2 cV-p 2 / CT + p 



2 



= C ^ dr 2 - dp 2 - ^ I In -— ^ 1 rffi . (33) 
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2. Flight velocity of galaxies and associated redshift 

The flight velocity of a galaxy at r = const can be obtained from eqs. ([19]) and (|29|) with 
a(t) =at. Alternatively, from eqs. ( [3T]) we immediately get 

olt p olt 

p = ct tanh — and v = — = c tanh — . (34) 
c r c 

1. For expansion coordinates, with a(i) = at the eqs. ( 1201) and (12~T]) yield 

- — and r pm = - — = - In 



> , — ■ em l , , 

whence 



Aem 

2. In explosion coordinates from eqs. (1221) and (1341) we get 

n 1/2 

1 + tanh(ar/c^ 



A " =e ar/c . (35) 



e Qr/c . (36) 



1 — tanh(ar/c) 

According to the first of the eqs. (1321) there is no additional gravitational redshift of the light 
on its way from p = p em to p = 0, i.e. \o/X p = 1, and therefore we obtain 

Ao_ = e ar/c ; (3?) 
''em 

exactly the same result as in expansion coordinates, eq. ( 135]) . It should be noted, however, 
that, following from eqs. (13~T|) . ([1]) and ( 133]) . in the two coordinate systems the distance from 
the point of emission as well as the time of emission differ from each other. 

3. The case of constant expansion velocity is particularly suited for demonstrating the 
compatibility of the two interpretations at issue, because in explosion coordinates the redshift 
is completely due to the Doppler effect. According to eqs. ( 1341 the spatial grid of expansion 
coordinates r is moving relative to that of explosion coordinates p at a speed that increases 
with the distance from the origin. We now consider a light ray emitted at r = r em and 
directed toward the origin. From the viewpoint of the system S exp i of explosion coordinates, 
observers at rest in the system S expa of expansion coordinates are moving relative to the 
emitting galaxy at r = r em at a velocity that increases with decreasing r. Therefore in S exp i 
they appear to observe a Doppler shift that continuously increases as r approaches zero. 
This way an observer in S exp i can readily understand why for observers in S expa the redshift 
appears to be accumulated on the way of light from emission to reception. 
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3. View of an external observer 



According to the theory of chaotic inflation [26ll27|. in the case k = (no spatial curvature) 
considered in this paper a Friedman-Lemaitre universe could be finite and embedded in 
an infinite super- universe. This so-called multiverse is filled with a "foam" of fluctuating 
quantum fields from which numerous universes of all kinds can emerge by inflation or have 
already done so. It is interesting to find out how the situation, so far considered from an 
internal observer, would be seen by an external observer located at some distance from our 
universe. 

Considerably simplifying the above model we consider a toy model of the universe in 
which a finite spherical section of a Friedman-Lemaitre universe is surrounded by a suffi- 
ciently large bubble of true vacuum, uninfiltrated by any externally generated gravitational 
fields. In contrast to the situation outside a collapsing star, due to Birkhoff 's theorem (time 
independence of all metric coefficients in the vacuum surrounding a spherically symmetric 



mass or energy distribution, see 28] or e.g. 29|]) the metric of the external space is not of 



the Schwarzschild type but must be pseudo- Euclidean. The reason is that it has been of this 
kind before the universe emerged by a creation out of nothing, a process tolerated by general 
relativity when the big bang (or, rather, an almost big bang) is preceded and triggered by 
an inflation field (dark energy, inflaton). Therefore, in polar spatial coordinates p, and <p 
the corresponding line element is 

ds 2 = c 2 dr 2 - dp 2 - p 2 dQ. (38) 

The case a(t) = at represents a fairly good approximation to the state of our universe from 
the end of the period of matter dominance until now, i.e. for the time interval to/2 < t < t 
where to is the present age of the universe. In this time interval the pressure can be neglected 
and matter can be treated as pressureless dust. The only forces acting on the matter elements 
are gravitational forces (including the action of dark energy), and we can therefore assume 
that in comoving coordinates the boundary of the universe (which must lie well beyond 
our horizon) is at rest, r = R = const, and is constituted by matter elements moving at the 
velocity (J34l) with r = R. 

Since according to eqs. ( |33i) and ( 1381) the external and internal coordinates are both not 
only Gaussian normal but also employ the same radial metric^ it is rather obvious that the 
2 For a smooth connection with the external coordinates, as internal coordinates also standard coordinates 
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velocity of the boundary seen by an external observer is the same as that seen by an internal 

observer. A more physical proof of this result is the following. We consider the propagation 

of light emitted with frequency vq at the origin and directed toward an observer outside the 

universe, at rest at p = p b s = const. According to the eqs. ( 133]) and ( )38|) . on its whole way 

inside and outside the universe the propagation of light is described by the equation p(r) = c 

or p = c(t—t ) respectively. In consequence the time between the arrival of successive wave 

crests at the place of the observer in "outer space" equals the time between their emission 

at the origin whence v obs = u Q . An observer comoving with the boundary of the universe sees 

the origin receding from him at the velocity v given by the second of eqs. ( 134|) . He therefore 

observes a redshift of the light from the origin, and according to eq. ( 12 2 j) he measures the 

frequency / 1 — v/c\ ^ 2 

\1 + v/c j 

The light leaves him with the same frequency as it had on its arrival. Since the boundary 
moves toward the external observer, the latter will observe a blueshift of the light from the 
boundary and measure the frequency 

(l + u/cV 12 (l-v/cV 12 (l+u/c\ 1/2 

where u is the velocity of the boundary in his (external) coordinate frame. As we said above 
we must have v b s = vo, and with this the last equation yields 

a R 



u 



c tanh (39) 



c 



where at last eq. ( 134)) and r = R was used. 



4- Numerical values of characteristic velocities 

1. In the comoving coordinates of the Robertson- Walker metric, according to Hubble's 
law the present recessional velocity of the boundary of the observable universe is v = H d 
where d — a(to) r bo w 3.5 ct is its metric distance from us and to the age of the universe. 
With Hq = 0.7 • 3.24 ■ 10~ 18 s _1 it becomes v rec = 3.45 c. A reasonable assumption about the 
(ds 2 = goo(p, t) c 2 dr 2 + g pp (p, r) dp 2 + p 2 dQ) could be envisaged. However, for a(t) ~ t, at some critical p 
a coordinate singularity occurs, and for all p > p cr u standard coordinates do not exist. This provides 
another example that compatibility with the requirements of general relativity may prevent the existence 
of coordinates with specified properties. 
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outer boundary R of the universe is that at present it is two times as far away from us as 
the present boundary of the observable universe, i.e. d ~ 7 ct . (This way inhomogeneities 
propagating from "outer space" into the universe cannot have spoiled the homogeneity and 
isotropy inside the observable universe, see [if]], page 231). The corresponding recessional 
velocity is 

v rec = Q.9c. (40) 

2. In explosion coordinates the velocity at which the present boundary of the observable 
universe moves away from us is given by the second of the eqs. (1341) . Following from a(t) =at 
we have a r bo /c = a(t ) r bo /(ct ) = d Q /(ct Q ) =3.5, and with this we obtain v bo = 0.998 c. At 
the outer boundary of the universe we have a R/c = d/(cto) = 7, and according to eq. ( 139]) 
the velocity at which it moves is 

u — v — 0.999998 c . (41) 



B. Universe with constant expansion rate 

We consider a second case in which a full analytic solution for explosion coordinates 
can be obtained, namely a(t) = a e Ht with constant expansion rate H. In this case which 
describes an inflationary expansion, in addition to the Doppler shift there is a gravitational 
blueshift of the light from far away galaxies. 



1. Derivation of explosion coordinates 

In the present case, instead of using eqs. (fT4"|) and (JIB]) it is easier to go back to eq. (J7|), 

t pp + Ht p = -H/c 2 . 

After solving the homogeneous equation, the solution of the inhomogeneous equation can 
be obtained by variation of constants and is given by 

1 Hp 

t(p, T) = t + — In cos w with w = . (42) 

H c 

An integration function was chosen such as to satisfy the second of the conditions ()6]). With 
this result and a(t) = ae Ht from eq. ( fTBl we get 

ce~ Hr 

r(p, r) = — — tan to. (43) 
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Inversion of the functions r(p, r) and t(p,r) given in eqs. fj42|) - (|43|) yields 



c 



aHr Tj,\ 1 , / aH 2 r 2 



" = n arcsin l^r eH, J • T = ( - 2E 'H 1 - V- H ■ (44) 

Inserting a(t) = a e Ht and eqs. fH2|) -f l43|) in eq. (TXTT) and the fourth of the eqs. (T4J) we obtain 

c 2 

^oo = cos 2 w , cfo = -— sin 2 it? . (45) 

From the first equation it follows that eq. ffl2|) is satisfied as it should according to the 
general theory. 



2. Flight velocity of the cosmic substrate and associated redshift 
From the eqs. (fT9j) . (|42p and a(t) = ae ffi we obtain the flight velocity 

v = c sinu> . (46) 



1. In the present case eqs. (1201) and (I21J) become 

An ^ .., c /■*<) ,. c 



A- = e H (*°-^) with r em = - (° e~ Ht dt = — (e~ Ht ™ - e~ Ht °) 



and in combination yield 

^em \ C ) 

From this, with r em — > r, t em — > t and a H r/c = e~ Hr tan to according to eq. (H3|) and using 



(47) 



eq. (]42]h for expansion coordinates we obtain the result 

Aq = 1 
A em 1 — sin w 

2. In the system of explosion coordinates, from eqs. (l22i) and ()46l) we obtain 

A p 1 + sin w 



X em cos w 
and from eq. ( 1241) and the first of the eqs. f j45l) we get 



(48) 



A /cos 2 I (4g) 



A p V cos 2 w cos w 



Combing the results (148]) and (|49l) we obtain the same result (|47|) as in the system of 
expansion coordinates. 
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3. Volume expansion rate 

According to eq. (jlSjl the volume of a spherical shell of small thickness Ap is given by 

c 2 

AV(p, t) = \.ix —- Ap sm 2 w . 
H l 

Because it is time independent, the expansion rate E = dAV/(AV y/goodr) vanishes every- 
where. This means that in the present case, exponential inflationary evolution, the cosmic 
substrate is moving relative to a completely invariable space. 



4- Acceleration of the cosmic substrate and comparison with ordinary explosions 



From the first of the eqs. (1441) and from eq. (l4*6l) with 



dp 
dt 



Hp 

c tan = c tan w 



c 

for simplicity expressed in terms of the time t we obtain 

dp 



... dw 
v(t) = c cos w — — 
y ' dt 



H cosw 

dt 



c H sin w = Hv 



From this there follows an exponentially growing acceleration of the cosmic substrate along 
its radial trajectories. Accordingly the inflation preceding an almost big bang describes the 
actual explosion whereas the subsequent dynamical processes are basically consequences of 
it. (In the case of a big bang without inflation the dynamics of explosion is compressed into 
a singular instant.) Like in ordinary explosions the accelerated substance - here dark energy, 
there e.g. a chemical explosive like dynamite or gunpowder - simultaneously is the blasting 
agent driving the explosion, and for an external observer it moves outward filling a previously 
empty space. What differs from ordinary explosions is that there are no shattered fragments 
of a casing, and there is no precursive shock wave. Instead there will be a precursive front 
of the inflation field, a weak discontinuity, which propagates at the speed of light. 

C. Hubble Law 

In expansion coordinates the distance from the origin of a galaxy or an element of the 
cosmic substrate is d = a{t) r. From this follows the Hubble law 

d(t) = Hd(t) with H = . (50) 
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In the explosion coordinates belonging to the case of constant expansion velocity, a(t) = at, 
according to the second of the eqs. (134]) we can write the flight velocity of galaxies etc. in 



the forms 1 

v = H expl p with H exp i = -. (51) 

T 

Formally, this is identical with the Hubble law 

d{t) = Hd(t) with H = - 

t 

obtained in expansion coordinates. Physically the two laws are different in that the under- 
lying lengths as well as proper times are measured differently. 

Let us now consider the case of constant expansion rate H. In expansion coordinates 
eq. ( 1501) with H = const holds. In explosion coordinates according to eq. (1461) we have 



v = H explP with H expl = y^ = #h __£- + ... j , (52) 

where at last H exp i was expanded with respect to Hp/c. In contrast to the result for ex- 
pansion coordinates, H exp i is weakly space dependent. This is not in contradiction to the 
result obtained for expansion coordinates or from observations although a space dependent 
Hubble parameter appears unusual. The reason is that measurements concerning far away 
objects involve the application of a theory related to the specific coordinate system in use, 
and this theory must appropriately be adjusted in the transition to explosion coordinates. 
(For example, the luminosity distance must be redefined.) 

III. APPLICATIONS AND EXTENSIONS 

1. In more general cases than the ones considered in the last section g 0Q will depend 
on p and r. Therefore, in addition to the Doppler shift a time dependent gravitational 
redshift of light will occur. While the Doppler shift can still be calculated from the Doppler 
formula ( l22|) . no generally valid formula for the shift effect of time dependent gravitation is 
known (at least to the author of this paper). However, for all cases covered by eq. ([2]) the 
latter can be derived by making use of the equivalence of the two different conceptions of 
galaxy recession. Combining eqs. (120]) . ( 121]) and (123]) . with t em — > t as in eq. (123]) we obtain 

A ° - A ° Aem - a{to) with r-r -cP — (53) 

Willi / — I em — C I j—r . K'-'d) 



\ P Km \ P a(r) (l+^/l-a 2 (t)/a 2 (r)) "" Jt em o(f) ' 
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2. Objects that are small in relation to cosmic distances like atoms, the solar system 
or even galaxies do not participate in the spatial expansion of the universe associated with 
usual theory. In expansion coordinates this phenomenon is not easily comprehensible. A 
first important proof of it is implicitly contained in one of A. Einstein's last papers [30]. It 
is shown there that a star in static equilibrium can smoothly be embedded in an expanding 
universe which means that the radius of the star does not expand. From this it is often 
derived that cosmic expansion is restricted to distances of cosmic scale. (According to 32 1 
"only distances between clusters of galaxies and greater distances are subject to the expan- 
sion".) What cannot be understood on this basis is, why, on the other hand, light rays of 
much smaller spatial extension are still subject to expansion. (This applies in particular to 
the incoherent light emitted by far away galaxies etc.. However, in the usual derivation of 
the cosmic redshift only monochromatic wave trains of infinite extension are considered.) 
This different behavior is especially difficult to understand in expansion coordinates 31], 
but it is very easily understood in explosion coordinates: Completely independent of the 
radial extension of a radially directed light ray the frequency of it is redshifted due to the 
Doppler effect and the gravitational field. 

3. The transformation from expansion to explosion coordinates can only in special cases 
be analytically expressed and will in general involve numerical calculations. A formulation 
of the general relativistic equations for the dynamic evolution of the universe in explosion 
coordinates may offer an easier approach to solutions, at least in special cases. 

4. The introduction of explosion coordinates in a universe with negative spatial curvature 
(k = —1) is certainly feasible in a similar way as in the case k = considered in this paper. 
In the case of positive spatial curvature the situation is different. Nevertheless it may be 
possible to successfully impose the condition g pp = also in this case, although this would 
appear somewhat artificial. 



IV. CONCLUSION 

It was shown in this paper that explaining the recession of galaxies (or other manifes- 
tations of matter at earlier stages) by an expansion of space or as an explosion and after- 
explosionlike motion relative to space is equivalent. Both interpretations are relative in that 
their validity is restricted to specific coordinate systems. The transition between them can 
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be performed by simple transformations. 

Certainly, the usual approach using expansion coordinates has the major merits. Most 
important is its simplicity provided by the fact that many physical properties are related 
to one simple time dependent parameter, the scale factor a(t). Furthermore, the coordinate 
time is equal to the proper time and is thus valid for the whole universe. For people unfa- 
miliar with the field the concept of an expanding space may occasionally provide difficulties, 
but this is by far surpassed by the afore-mentioned advantages. 

Our approach by explosion coordinates aimed at an alternative interpretation and conse- 
quently is based on a reformulation of known results. Nevertheless it has its merits as well. 
For one, it better fits the view of an external observer. Also, the interpretation of the big 
bang or inflation as giant explosions, and the restriction of recessional velocities to values 
below the velocity of light are more intuitive. Furthermore there exist problems which are 
easier to handle in explosion coordinates. Finally, the close affinity to Newtonian cosmology 
may provide advantages in some cases. 

The two interpretations should by no means be confounded. As well as the Doppler 
effect must not be used for explaining the redshift in expansion coordinates, their corre- 
sponding distances and recessional velocities must not be employed when Doppler effect and 
gravitational blueshift are evaluated in explosion coordinates. 
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